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We study the effect of the charm quark mass in the CTEQ global analysis of parton distribution 
functions (PDFs) of the proton. Constraints on the MS mass of the charm quark are examined at 
the next-to-next-to-leading order (NNLO) accuracy in the S-ACOT-x heavy-quark factorization 

rf^ \ scheme. The value of the charm quark mass from the hadronic scattering data in the CTIO NNLO 

fit, including semiinclusive charm production in DIS at HERA collider, is found to agree with the 

f^ ' world average value. Various approaches for constraining rUc in the global analysis and impact on 

^SJ I LHC cross sections are reviewed. 
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I. INTRODUCTION 

r^ , Quark masses are free parameters of the QCD Lagrangian that parametrize explicit breaking of the chiral symmetry. 

^-t' For quarks heavier than 1 GeV, quark masses arise as independent hard scales mg in perturbative QCD calculations 

Oj! for particle cross sections. Quarks arc not observed freely because of color confinement, hence their properties and 

QJ • masses are established indirectly by comparing theory calculations against experimental data on hadronic reactions. 

1-^ ' In a global analysis of parton distribution functions (PDFs) of the proton, the method by which the heavy-quark 

masses are included in experiments at energies comparable to mg has non-negligible impact on the extracted PDFs 

[l|. Variations in the PDFs associated with the treatment of heavy quarks have phenomenological consequences for 

^ , electroweak precision measurements at the Large Hadron Collider [^ . 

^^ ' In this paper we explore constraints on the charm quark mass from the global hadronic data in the CTEQ NNLO 

0^ " PDF analysis. The study is motivated by publications of combined cross sections on inclusive deep-inclusive scattering 

'^ I (DIS) and semi- inclusive DIS charm production at the ep collider HERA [3, |j]- Among all experimental data sets 

included in the global fit, the DIS experiments have the best potential to constrain the charm mass. On the theory side, 

■^ I PQCD calculations for neutral-current deep inelastic scattering have been extended to the 2-loop level in the QCD 



X 



\r I t^ycu calculations tor neutral-current aeep melastic scattermg have been extenaea to the z-loop level m the UiJ 
^^ ■ coupling strength a^ both for massless 0-01 and massive |8l-[l0| quarks, while massless [ii|,[i3l and some massive [li. 
^_2 . [201 coefficient functions were also obtained at the 3-loop level. With such accuracy, it becomes possible to determine 
the charm quark mass and its uncertainty from the DIS data. 

The current world average of the running charm mass defined at two-loop QCD in the MS renormalization scheme is 
'mcirric) ~ 1.275 ±0.025 GeV [21[. This value is derived primarily from measurements in timelike scattering processes 
and lattice simulations. The analysis of the DIS data from HERA in principle allows us to extract mc{mc) from a 
C^ . spacelike scattering process. 

In many previous PDF analyses, the heavy-quark masses have been treated as effective parameters rather than 
fundamental constants. They were anticipated to deviate from the MS masses or even be fully independent. Besides 
entering the exact matrix elements, the heavy-quark masses control various approximations in DIS calculations, as 
will be reviewed below. The approximations affect extraction of the mass values from the hadronic data, but their 
effect is of a higher order in the QCD coupling strength according to the QCD factorization theorem. Thus, although 
the bias due to the approximations can be important at low orders, it is expected to subside as more loops are included 
in PQCD calculations. At high enough order of as, such as the NNLO, a comparison to the world-average MS quark 
mass becomes feasible. 

There exist several theoretical approaches, or heavy-quark schemes, for computations involving massive quarks 
[23 - I28J . The extracted rric value depends on the heavy-quark scheme as well as the order of the PQCD calculation. 
In this work, we adopt a two-loop implementation [23| of the general-mass scheme S-ACOT-x [23, [30, ^M employed 
in the CTIO NNLO analysis [32|- We discuss physics assumptions behind the S-ACOT-x scheme and differences 
compared to the other schemes used in the recent extractions of rric [j, [SJ, [SJI • Within this scheme, we find the range 
of input mc{mc) values providing the best description of the CTIO fitted data and compare the accuracy of extraction 
to the uncertainty on the world-average value. Finally, we analyze the impact of the uncertainty in the charm mass 
on benchmark LHC predictions. 



II. CHARM MASS DEFINITIONS AND NNLO PREDICTIONS FOR DIS 

A. Overview 

Quantum field theory operates with two common definitions of the quark mass, the pole mass and the MS mass. The 
pole mass is defined as the position of the pole in the renormalized quark propagator. The pole mass is infrared-safe, 
gauge- invariant, and is derived in the on-shell renormalization scheme. As the charm mass is close to the border of 
the nonperturbative region, the accuracy of the pole mass determination for charm is limited by significant radiative 
contributions associated with renormalons (35l - l37l | . Because of large perturbative coefficients arising even at three or 
four loops in the QCD coupling a^, the pole rric value cannot be determined to better than a few hundred MeV. 

The MS mass mc{^), on the other hand, is the renormalized quark mass in the modified- minimal-subtraction 
scheme, defined as a short-distance mass that is not affected by nonperturbative ambiguities. It is evaluated at a 
momentum scale ^ typical for the hard process, often taken to be the mass rric itself. Precise determinations of 
mc{mc) achieve a smaller uncertainty of order 30 MeV or less. The MS mass starts to differ from t he p ole mass 
beginning at order 0{as)- The conversion between the MS mass to the pole mass may be required J38l - l4l| . because 
the experimental mass definition is often closer to the on-shell mass. 

In parton-level diagrams for deep-inelastic scattering, external massive quarks may arise both as the final and initial 
states. For the quarks that are heavier than the proton, some factorization schemes introduce an effective PDF to 
describe their quasi-collinear production at high energy. The heavy-quark PDF can contribute to the hadronic cross 
section through a convolution with a hard-scattering matrix element with the heavy quark(s) in the initial state, 
also called a "fiavor-excitation" matrix element. In contrast, the "flavor-creation" matrix elements include only light 
quarks and gluons in the initial state, while the heavy quarks are only in the final state. The "flavor-excitation" 
contributions commonly arise in the variable flavor number (VFN) schemes, such as the general-mass VFN (GM 
VFN) scheme. The alternative fixed-fiavor number (FFN) scheme does not introduce a heavy-quark PDF operates 
with "flavor-creation" terms only. 

For this study we employ the S-ACOT-x general-mass scheme [23, |30, E^] implemented to the 2-loop (NNLO) 
accuracy |29[ . This scheme includes exact massive flavor-creation contributions that dominate at low boson virtualities 
Q, as well as the approximate flavor-excitation terms that are important at high Q. Thus, the S-ACOT-x scheme 
reduces to the FFN scheme at Q^ « m^ and to the zero-mass VFN scheme at Q^ » m^. 

In a comprehensive factorization scheme such as GM-VFN, the charm mass plays several roles. First, the exact 
charm mass enters Feynman diagrams for charm particle creation in the final state, such as 7*3 — > cc in NC DIS. 
Second, auxiliary scales are introduced that are of order of the fundamental charm mass (either the pole mass or MS 
mass), but need not to coincide with it. 

One such scale sets the energy for switching from the 3-flavor to 4-flavor evolution in the running asip-)-, which is 
utilized by both FFN and VFN computations. A similar switching scale from 3-flavor evolution to 4-flavor evolution 
arises in the PDFs f g/pjx , /i). The charm mass also deflnes characteristic energy scales in the flavor-excitation 
contributions, cf. Sec. IIIBI Finally, there can be auxiliary scales associated with the final-state quark fragmentation 
into hadrons, present both in the FFN and GM-VFN schemes. The dependence on these auxiliary scales is reduced 
with each successive order of perturbation theory. 

When the input mass is varied in the global fit, the response of the DIS cross sections reflects coordinated variations 
of all such scales. An important question arises when interpreting the outcome of the fit: which mass parameter 
controls the agreement with the data, the exact charm mass or the approximate mass in the auxiliary scales? 

We have found that the global fit is sensitive to the exact charm mass, despite the introduced approximations. In 
one exercise, we have varied the input charm mass in the exact DIS coefficient functions for fiavor-creation processes, 
while keeping it fixed in the above auxiliary mass scales. In a complementary exercise, we varied the input mass in all 
auxiliary scales, while keeping it fixed in the exact fiavor-creation coefficient functions. In both cases, we examined 
the agreement with the data as a function of the varied mass parameter. We followed the fitting procedure outlined 
in the next section and assumed fixed PDF parametrizations for the best- fit mc{mc) found in the main analysis. 

The dependence of the figure-of- merit function x^ on the varied mc(rac) in these exercises is shown for all fitted 
experiments, the combined HERA inclusive DIS data [3|, and the combined HERA semi-inclusive charm production 
data [1] in the upper left, upper right, and lower panels of Fig. [TJ To better visualize the comparison, x^ is divided 
by the number Npt of data points for each data set. The solid blue line and dashed magenta line are for x^ /Npt for 
the varied mass parameter in the exact DIS coefficient functions and in the auxiliary mass scales, respectively. 

With all PDF parameters fixed, variations in x^ are more pronounced in the mass scan of the first type, when the 
charm mass is varied only in the exact DIS coefficient functions. In this case, x^ in both inclusive and semiinclusive 
DIS shows a pronounced minimum as a function of mc{fnc)- 

In the second case, when mc{mc) is varied only in the auxiliary scales, the x^ dependence is flatter and has a 
shallow minimum at most. This exercise indicates that both inclusive and semiinclusive DIS cross sections are more 



sensitive to the exact rric mass in the flavor-creation cocfBcicnt functions than to the auxiliary scales. The detailed 
constraints on rric are determined by the interplay of these two dependences as well as by variations in the PDFs and 
other inputs of the global fit. 
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FIG. 1: Dependence of y^ /N^t as a function of mc{mc) in the exact flavor-creation coefficient functions (solid blue lines) and 
auxiliary energy scales listed in the text (dashed magenta lines). 



B. Details of implementation 



Our calculation proceeds by taking the MS quark masses as the input for the whole calculation. The transition 
from the 3-flavor to 4-flavor evolution in ag and PDFs is taken to occur at the scale equal to this input mass.* 

The massive 2-loop coefficient functions for neutral-current DIS with explicit creation of cc pairs [9| and the operator 
matrix elements AJ^^ [2j] that we use require the pole mass as their input. For these parts, the MS mass is converted 
to the pole mass according to the 2-loop perturbative relation in Eq. (17) of [4J]. In charged-current DIS, for which 
only 1-loop expressions are included, we can use either the MS mass or pole mass and find almost no sensitivity to 
the choice. 

The general-mass schemes |23l428j differ primarily in the form of approximation for flavor-excitation coefficient 
functions (with incoming heavy quarks) at Q comparable to rric- When rric is negligible, the coefficient functions 
in all GM-VFN schemes reduce to the unique zero-mass expressions, but, near the threshold, they may differ by 
powerlike contributions (TTi^/Q^)^ with p > in the approximate flavor-excitation terms. The powerlike corrections 
are suppressed in the full result by an additional power of as, i.e., they are of order a^ in an NNLO (aj.) calculation. 



The evolution of as and PDFs is carried out using the HOPPET computer code [43 
flavors occur at the MS masses. 



, configured so that transitions from A''^^ to Nf^i 



Numerically, the powerlike contributions can explain the bulk of the differences between the GM VFN schemes, 
whether in the extraction of rric or other computations. 

In the S-ACOT-x scheme, the form of the powerlike contributions is selected based on the general consideration of 
energy- momentum conservation j3ll |42| | . The requirement that the flavor-excitation contributions obey the conser- 
vation of energy suppresses them at energies close to the massive quark production threshold and results in better 
description of the DIS data. The flavor-excitation coefficient functions are constructed from the respective zero-mass 
expressions evaluated as a fraction of the rescaling momentum fraction x(iVf^) that depends on the total mass Alf 
of heavy quarks in the final state. The general behavior of x is determined from the condition of the threshold 
suppression, while its detailed form can be varied to estimate the associated higher-order uncertainty in the extracted 
rric- This approach is readily demonstrated [23| to be compliant with the QCD factorization theorem for DIS cross 
sections to all orders [43, which is one of the advantages of using the S-ACOT-x scheme. 

In the default implementation, the rescaling variable in charm production is given by x = 2;(1 -I- Arri^/Q^), corre- 
sponding to Mf = 2mc for the lightest final state {cc). ^ The rescaling ratio x/^ is thus independent of x. On the 
general grounds, it may be expected that the threshold suppression is less pronounced at W^ = Q'^{l/x — 1) — > 00 for 
fixed Q, corresponding to x — >■ 0. In this limit, it may be desirable to reduce or even eliminate the rescaling altogether, 
as quasi-coUincar production of heavy quarks becomes more feasible. 

To allow for this possibility, a generalized rescaling variable C, can be implicitely defined by [42| 

where A is a positive parameter, typically < A < 1. The S-ACOT-x scheme is reproduced with A = 0, and the 
rescaling is fully turned off for A ^ 1. For A 7^ 0, the mass-threshold constraints are enforced at small W (large x), 
but the standard x variable is recovered at large W (small x) in a controllable way. 

The CTEQ global fit is sensitive to variations in A. The changes in the preferred value of rric that we observe 
provide an estimate of the uncertainty due to the powerlike corrections. 

III. CHARM MASS SCANS 

Using the theoretical setup reviewed in the previous section, a scan of the log-likelihood function x^ over the input 
MS charm mass in the range 1 < mcimc) < 1.4 GeV was carried out by utilizing the setup of the CTIO NNLO PDF 
analysis |32| . Besides the combined HERA data on inclusive DIS and semiinclusivc charm production, we include 
experimental data from DIS measurements by BCDMS ]M,kA: NMC Ei, CDHSW 0, and CCFR (MIU; NuTeV 
and CCFR dimuon production l53,J53|; ^20 measurements at HERA [5J| that are not included in the combined set; 
fixed-target Drell-Yan process J55l - [57| : vector boson and inclusive jet production at the Tevatron [58l - l65| . We also 
include inclusive jet production at the LHC [63, [631, which slightly reduces the uncertainty in the gluon PDF. 

In the several series of fits that we have carried out, we varied the assumptions and inputs to test their impact on 
the optimal charm mass. The results that we will present are common to all these fits, hence are robust against the 
variations in the assumptions. 

Depending on the candidate fit, the QCD coupling strength was taken to be either the world average, as{Mz) = 
0.118 [2l[, or a lower as{Mz) = 0.115, which is preferred by the CTIO NNLO analysis when as{Mz) is allowed 
to vary. The factorization/renormalization scale /i in DIS was set equal to the momentum transfer Q. To test the 
sensitivity to the PDF parametrization form, the initial scale Qq at which the input PDFs are provided was either set 
to be Qo = 1 GcV independently of rUc or varied in the scan together with rric as Qq = rric — 0.005 GeV. Several forms 
of the gluon PDF parametrization were considered, since DIS charm production is sensitive to the gluon PDF g{x, Q). 
At the initial scale Qoi we either constrained g{x, Qq) to be positive at all x or allowed it to be negative at small x, 
provided the negative gluon did not lead to unphysical predictions. In the latter case of the negative gluon, the fit 
included the HI data on the longitudinal structure function Fl{x, Q) |68| and an additional theoretical constraint to 
ensure positivity of Fl{x, Q) at x > 10^^. 



t Starting from 0{a'^), contributions with up to four massive quarks in the final state can appear. In such terms, we still use x 
a; (1 -I- 4m^/Q^), given their smallness in the total result 1291 . 
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FIG. 2: Global ^ of the S-ACOT-x NNLO fit as a function of the MS charm mass. Lines with left/right arrows indicate 90% 
C.L. intervals obtained with different tolerance criteria. 

A. Results of the fit and error analysis 

A typical plot of the global x^ in the above mass scans as a function of mj^m^ is shown in Fig. [51 In this plot, the 
PDF parameters were refitted for every 1x1^ value so as to minimize the global x^ compared to the CTIO NNLO fit. 
The default rescaling variable x was used, corresponding to the rescaling parameter A = 0. 

In the figure, the preferred value of the MS charm mass, corresponding to the minimum of x^, is mc(rn^ = 1.12 
GeV. The optimal charm mass is below the world-average mc{mc) = 1.275 ± 0.025 GeV, and this general trend is 
observed in all fits. 

From the mass dependence of x^ in the figure, we can compute the uncertainty on this value due to the PDF 
parameters according to several criteria available in literature. To define the 90% confidence level (90% C.L.) on 
"T-c(wc), we compared 3 different criteria for defining the PDF uncertainty: 1) the uniform x^ tolerance, in which 
one assigns a 90% C.L. to a Ax^ < 100 variation as in the CTEQ6 PDF analysis [11]; 2) the CTlO-like criterion, 
which supplements the uniform x^ tolerance condition by additional x^ penalties to prevent strong disagreements 
with individual experiments on average j70| : 3) and the MSTW-like criterion, which does not introduce the uniform 
tolerance, but requires the x^ value for every individual experiment to lie within the specified confidence interval. In 
the latter two methods, deviations of the PDF parameters from their best-fit combination are additionally constrained 
so as not to trigger a strong disagreement with one of the fitted experiments. This condition effectively reduces the 
PDF uncertainty in methods 2 and 3 compared to the uniform tolerance criterion (method 1), see further details in 

The PDF uncertainties on mc(mc) obtained with the three definitions are reported in Table ID The uncertainty 
according to the uniform tolerance definition is larger than with the other two definitions, as anticipated. 



PDF uncertainty Snic [GeV] (90% C.L.) 


Ax^ < 100 


5m. =t«:i° 


CTlO-like 


5mc =±ri^ 


MSTW-like 


x^ +0.12 

OTUc =_o.l8 



TABLE I: The PDF uncertainty on the optimal value of vac extracted from the fit at the 90% C.L. by using three different 
criteria. 



Among all experiments included in the mass scan, the tightest constraints on rric are imposed by the inclusive DIS 
HERA data [31 and HERA charm production [J]. Weaker constraints also arise from the other DIS experiments. 
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FIG. 3: Individual x^ contributions of the combined HERA-1 inclusive DIS data, NMC F2 data, and combined HERA charm 
quark production data as a function of mdrac). 



notably the NMC measurement of -Ff (x, Q). 

Fig. 13] shows x^ as a function of rric in these three experiments in a representative mass scan. The inclusive DIS 
data from HERA (upper left inset) and NMC (upper right inset) broadly prefer an rric in the range 1.1-1.3 GeV, while 
the charm HERA data (upper right inset) prefers lower rric values of order 1 GeV. The NMC shown in the lower inset 
prefers a higher nic of order 1.2 GeV in accord with the inclusive HERA data, but with flatter y^ dependence. While 
in the shown mass scan the inclusive DIS data prefer a higher nic value than in DIS charm production, the minimum 
of x^ in inclusive DIS may shift to lower values of about 1.05-1.1 GeV depending on the parametrization forms of the 
PDFs and other fit assumptions. 

We note that the published inclusive reduced cross sections on HERA charm production were derived from charm 
differential distributions by applying a significant acceptance correction computed with the program HVQDIS [711] in 
the fixed-flavor number scheme. Due to the mismatch in the schemes, the mc{mc) preferred by the HERA charm data 



may be biased when determined in the S-ACOT-^ general-mass scheme. Nevertheless, we see from Fig. |3] that the 
mc(r?T-c) values that are separately preferred by the inclusive DIS and charm DIS data sets are generally compatible. 
In another cross check, we used separate, rather than combined, data sets for HERA charm production [72 - l79l |. The 
mdi^c) central value and its uncertainty from this fit did not change significantly compared to the fit to the combined 
HERA charm data set. 



Theoretical systematic uncertainty 


DIS scale 
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X definition 
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[Q/2, 2Q] 


[0.116, 0.120] 


[0, 0.2] 
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Smdnic) (GeV) 


+0.02 
-0.02 


+0.01 
-0.01 


+0.14 
-0 


+0.06 
-0 



TABLE II: Shifts of the optimal value of the charm mass inc{mc) obtained by varying theoretical inputs. 



Table |IT] summarizes systematic uncertainties on rric associated with theoretical inputs, including the factoriza- 
tion/renormalization scale, as{Mz), the A parameter in the rcscaling variable, and implementation of experimental 
correlated systematic errors. In this list, the last source of uncertainty arises from the existence of several prescriptions 
for including correlated systematic errors from the fitted experiments into the figure-of-merit function x^ [32, ISO] ■ De- 
pending on the prescription, the relative correlated errors published by the experiments can be interpreted as fractions 
of either central data values or theoretical values. This leads to numerical differences in absolute correlated errors, 
which may affect the outcomes of the fit, as discussed in the above references. As a default choice we normalize the 
correlated systematic errors to the theoretical predictions and estimate the uncertainty by changing the normalization 
to the experimental data for the DIS processes. These methods are designated as "extended T" and "D" methods in 
Ref . [32| . The uncertainties on rric associated with the parameter A in the rcscaling variable and definition of x^ a-re 
dominant. 

To further illustrate the impact of these two uncertainties, we show in Fig. |4] a contour plot of the dependence of 
the extracted mc{mc) on the A parameter and definition of x^. The values of A adopted in the mass scan are indicated 
on the horizontal axis. The triangle and star symbols indicate the best-fit values of mc{mc) obtained with the default 
and alternative prescriptions for implementing the correlated systematic errors in x^- The y^ minima are reached 
when A sa in both cases, and the best- fit mc{mc) values are 1.12 and 1.18 GeV, respectively. 

The solid red and black dashed contours correspond to the 68% and 90% C.L. PDF uncertainties for the two x^ 
prescriptions. -I- They overlap with the world-average interval mdmc) = 1.275 ±0.025 GeV [2l|, indicated in the figure 
by the horizontal blue band. We see that, by allowing a moderate increase of x^ compared to the best-fit value, and 
possibly a small increase in A, one can reconcile the optimal TOc(fTic) from the fit with the world-average value of 
mc{mc). 
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FIG. 4: Preferred regions for mc{mc) vs. the rescaling parameter A. The best-fit values and confidence intervals are shown for 
two alternative methods for implementation of correlated systematic errors. 



■f The CTIO or MSTW-like tolerance criteria lead to about the same boundaries. 
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FIG. 5: Relative changes in select PDFs fa/pi^, Q) at Q = 2 GeV obtained in a series of PDF fits with mc{mc) ranging from 
f to 1.36 GeV, plotted as ratios to the respective PDFs for mc{mc) = 1 GeV. Default rescaling (A = 0) in DIS coefficient 
functions is assumed. Darker colors correspond to larger mc{mc) values. 



B. Implications for PDFs and collider observables 



As 771c (ttIc) changes in the mass sean, parametrizations of the PDFs are adjusted so as to maintain agreement with 
the data. Representative best- fit PDFs from a charm mass scan are plotted in Fig. [5] at Q = 2 GeV. In this example, 
we vary the charm mass in the interval 1.0 < rac{rac) < 1.36 GeV that is about the same as the CTlO-likc PDF 
uncertainty determined in the previous section. For each mc{mc) selected from this range, we refit the PDFs, while 
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FIG. 6: Plot of NNLO cross sections for W , Z , Higgs boson production through gluon fusion, and top quark pair production 
at the LHC (8 TeV) for charm quark mass mdmc) ranging from 1 to 1.36 GeV and A = {0, 0.02, 0.05, 0.1, 0.15, 0.2}. Darker 
color corresponds to larger mass values. The empty diamonds represent the cross sections evaluated by using mc{mc) = 1.28 
GeV. The triangle and ellipse indicate the central predictions and 90% C.L. interval based on CTIONNLO predictions. 



keeping the other theoretical inputs at their default values, in particular assuming aa{Mz) — 0.118 and A = in the 
rescaling variable. Darker color points in Fig. [5] correspond to larger mass values. 

In general, as the charm mass is increased, both the charm PDF (upper left subfigure) and charm contributions 
to DIS cross sections are suppressed. Consequently, the gluon PDF (upper right subfigure) is enhanced in the 
intermediate x region, 10~^ ~ 10~^ so as to partly compensate for this reduction. This is accompanied by moderate 
enhancements in the up and down quark PDFs (lovifer subfigures) at x sa 10^^ — 0.5 and slight suppression of the 
same PDFs at a; sa 10"-^. 

PDF variations of such magnitude may have impact on collider observables (3|- For example, in Fig. [5] and [7] we 
show the dependence of NNLO total cross sections for W, Z, Higgs boson production through gluon fusion, and top 
quark pair production at the LHC at \S = 8 and 14 TeV. The NNLO cross sections for W and Z production are 
computed using FEWZ2.1 [8l|,|8^. The NNLO cross sections for Higgs and top quark pair production are obtained 
from iHixl.3 [sl and Top++2.0 [H [H with mn = 125 GeV, m* = 173.3 GeV, and the QCD scales set to the 
corresponding mass values. For each pair of total cross sections, we show the central CTIONNLO prediction and 
an ellipse corresponding to the 90% C.L. PDF uncertainty interval of the CTIONNLO set [331 ■ In the same figure, 
scattered points indicate the cross sections obtained with various combinations of mdrnc) and A in the intervals 
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FIG. 7: Same as in Fig. H but with \/S = 14 TeV. 



1 < mc{mc) < 1.36 GeV and < A < 0.2. Here the darker colors represent larger values of rric, as in Fig. [SJ The W , 
Z, and Higgs production cross sections increase with the charm mass by a few percent in the mass range considered. 
The trend is opposite for ti production. The changes are contained within the CTIONNLO PDF uncertainty ellipse, 
however, and do not modify the PDF-induced (anti-)correlations observed between the shown total cross sections. 

The shown uncertainty in the LHC cross sections due to rric is comparable to the experimental PDF uncertainty 
and in principle should be included independently from the latter. Let us outline one possibility for reducing the 
rric uncertainty. Instead of allowing mdmc) to vary in the whole interval 1 — 1.36 GeV of its PDF uncertainty, we 
could set it to be at the world-average mdmc) value of 1.275 GeV or in the la interval ±0.025 GeV around it. This 
would reduce the associated uncertainty in the PDFs and QCD observables. The corresponding predictions for the 
LHC W, Z, Higgs, and top quark pair production cross sections, obtained for mc{mc) = 1.28 GeV and five explored 
A values, are shown in Figs. [6] and [7] by open squares, with the size of the squrare increasing with the value of A. The 
spread in these predictions constitutes only a part of the full span covered by the scattered points for the interval 
1 < rric < 1.36 GeV. It reflects only the uncertainty due to the form of the rescaling variable, controlled by the A 
parameter. Theoretical predictions are better clustered in this case. 
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IV. CONCLUSIONS 

We explored the charm quark mass dependence in the CTEQ NNLO global PDF analysis that includes the recently 
published combined data on charm quark production at the ep collider HERA. This analysis, carried out in the 
S-ACOT-x heavy-quark scheme, renders an optimal MS charm mass that is compatible with the world-average value 
rricirnc) = 1.275 GeV. The best-fit value and its 90% C.L. PDF uncertainty (defined according to the CTIO criterion) 
are mc{mc) = 1.12lQ'i[7 GeV. Other significant sources of uncertainty exist and arc listed explicitly in Table HIl with 
the quadrature sum of their estimates given by Iq 02 GeV. 

In QCD predictions for massive-quark DIS, one draws a distinction between the MS charm mass, the fundamental 
parameter of the QCD Lagrangian, and auxiliary energy scales of order of the physical charm mass. The auxiliary 
scales that can contribute are specified by the factorization scheme. They can be associated with the evolution of the 
QCD coupling strength, evolution of PDFs, heavy-quark fragmentation, and powerlike contributions {rric/QY in DIS 
coefficient functions with initial-state charm quarks. 

We argue that the sensitivity to the auxiliary scales is reduced as the order of the PQCD calculation increases. In 
support of this argument, Fig. [T] demonstrates that the DIS data are sensitive mostly to the physical mass parameter 
in the exact DIS coefficient functions and less to the auxiliary mass scales in the other parts of the calculation. Thus, 
the hadronic cross sections in DIS and other processes at NNLO and beyond become increasingly suitable for the 
determination of the fundamental charm mass. 

The main findings of our fit are summarized by Fig. 21 showing the central values of the MS charm mass and 
contours of its PDF uncertainty as a function of a theoretical parameter A controlling the behavior of DIS coefficient 
functions with initial-state charm quarks. The central value of mc(mc) in the S-ACOT-x scheme tends to undershoot 
the world-average value according to the figure, but is compatible with the latter within the uncertainty. 

The uncertainties in the mcirric) determination of both experimental and theoretical origin were explored in Sec. lIIII 
If mc{mc) is varied as an independent parameter in the full range of order 1-1.4 GeV allowed by the NNLO PDF fit, 
it increases the uncertainty on the resulting PDFs, compared to a fixed vric- 

For comparison, the accuracy of the world-average mdmc), at about 0.025 GeV, is better by an order of magnitude 
than the NNLO fit uncertainty. By using a constant value of the mcirric) parameter, for example by setting it equal 
to its world-average value, one can suppress the corresponding uncertainty in the PDFs. This strategy is similar 
to implementing the QCD coupling dependence in the global analysis [86[, when using the world-average value of 
as{Mz) = 0.118 results in tighter constraints on the PDFs than in a fit with a free as{Mz)- When the input mdmc) 
value coincides with the world-average value, instead of being fitted, one suppresses the associated uncertainty in the 
PDFs and LHC cross sections, see sample calculations in Figs. [5] and [T] The residual theoretical uncertainty for a fixed 
mc{mc) then arises only from variations in the auxiliary scales and powerlikc contributions and can be suppressed by 
including higher orders in as- 
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